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Measurements  of  the  Effectiveness  of  Concave  Spherical 
Dimples  for  Enhancement  of  Hot-gas  Side  Heat  Transfer 


Edward  B.  Coy  and  Stephen  A.  Danczyk1 
Air  Force  Research  Laboratory,  Edwards  AFB,  CA,  93524 

Nomenclature 

d  =  Dim  pie  diameter 

h  =  Dimple  depth 

H=  Channel  height 

/  =  Fraction  of  surface  covered  by  dimples 

t  =  Dimple  spacing 

Nu  =  Nusselt  number 

Pc  =  Cham  her  pressure 


I.  Introduction 

In  a  n  ex  pander  cy  cle  roc  ket  engi  ne,  t  he  t  urbo-pump  i  s  dri  ven  by  a  g  aseous  p  ropellant  t  hat  has  been  heat  ed  i  n 
channels  incorporated  into  the  thrust  chamber  wall.  Th  e  power  output  of  the  turbine  is  limited  by  the  hot  gas  wall 
heat  transfer  c  oefficient  and  surface  area.  To  increase  the  chamber  pressure,  additional  pump  power  and  heat  a  re 
required.  Part  of  the  increase  can  come  from  an  increased  heat  transfer  coefficient  at  elevated  chamber  pressure,  but 
this  is  not  sufficient;  and  engine  designs  have  also  required  an  increase  in  surface  area,  typically  by  e  longating  the 
chamber,  or  a  supplemental  increase  in  heat  transfer  coefficient  by  one  of  the  methods  of  heat  transfer  enhancement. 
In  the  Orbital  Transfer  Vehicle1  program  several  styles  of  longitudinal  ribs  were  examined,  and  enhancements  of  40- 
50%  were  measured;  however,  there  are  other  approaches  which  should  also  be  considered  because  the  tips  of  ribs 
are  vulnerable  to  hot  spots  and  may  become  points  of  incipient  failure. 

1  Propulsi  on  Research  Engineer,  Aerophysics  Branch,  Space  and  Missiles  Propulsion  Division,  Propulsi  on 
Directorate,  10  East  Saturn  Blvd. 
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In  this  study,  we  examined  the  effectiveness  of  concave  spherical  dimples  (CSD).  T  his  method  has  been  the 
subject  of  several  previous  investigations  and  has  been  found  to  have  favorable  characteristics  with  respect  to  th  e 
efficiency  of  heat  transfer  enhancement  relative  to  the  amount  of  skin  friction  created.  CSD  may  also  be  less  prone 
to  the  formation  of  hot  spots  and  c  oncentrations  o  f  mechanical  stress.  The  main  objective  of  t  his  study  was  t  o 
develop  a  design  space  for  heat  transfer  enhancement  by  characterizing  the  heat  transfer  effectiveness  as  a  function 
of  non-dimensional  geometric  p  arameters.  Th  is  info  rmation  sho  uld  be  of  utility  to  th  rust  ch  amber  design  teams 
when  performing  trade  studies  involving  heat  transfer  enhancement  with  other  critical  factors  such  as  performance, 
weight  and  lifetime. 

CSD  patterns  can  be  characterized  with  two  non-dimensional  numbers:  the  depth-to-diameter  ratio,  h/d,  and  the 
fraction  of  the  surface  covered  by  dim  pies,  f  Ot  her  geometric  par  ameters  t  hat  may  be  si  gnificant  i  n  som  e 
applications  are  the  cha  nnel  height,  H,  a  nd  the  arrangement  of  th  e  r  ows  of  di  mples,  staggered  or  in-line.  For 
reference,  the  maximum  possible  value  of/is  0.91  and  occurs  when  the  edges  of  adjacent  dimples  are  tangent. 


Previous  works  on  heat  transfer  enhancement  using  CSD  are  reviewed  and  summarized  in  Li  gran  r  (2003). 
The  majority  of  previous  efforts  have  been  performed  using  nearly  ambient  air  in  flat  channels  with  H/d  in  the  range 
0.2-2.  Moon3  (2000)  examined  the  effect  of  channel  height  on  Nusselt  number  for  a  staggered  pattern  of  CSD  with 
h/d=0.19  and  /M). 57.  The  globally  averaged  heat  transfer  enhancement,  Nudimpis/Nusm0oth  was  approximately  2.1  and 
independent  of  channel  height  for  a  ra  nge  of  H/d  of  0.35-1.5.  M  ahmood4  (2002)  described  the  flow  patterns  and 
measured  local  Nusselt  number  distributions  over  the  surface  for  the  same  dimple  pattern  used  by  Moon.  The  most 
important  feat  ures  of  t  he  fl  ow  c  onsisted  of:  ( 1 )  a  shea  r  1  ayer  t  hat  separat  es  f  rom  the  u  pstream  edge  a  nd  t  hen 
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reattaches  on  the  downstream  edge,  (2)  vortices  that  shed  from  the  indentations  and  then  advect  over  the  flat  surface 
downstream  of  t  he  dimples,  and  (  3)  peri odic  unst  eadiness  as  fl  ow  i  s  ej  ected  and  t  hen  i n-mshes.  Local  Nus  selt 
numbers  were  in  the  range  of  1-3.2  times  those  for  a  smooth  plate.  Burgess5  (2003)  examined  the  effect  of  varying 
h/d  from  0. 1-0.3  for  a  fixed  value  of/=0.57  and  a  fixed  channel  height,  H/d=l.  The  enhancement  increased  with  h/d 
and  was  correlated  with  the  following  equation: 


NUfomp  le 

Nil  smooth 


=  1  +  6.193 


h_ 

\dj 


(l) 


Burgess  also  examined  the  e  ffect  of  turbulence  intensity  in  the  channel  and  found  that  it  had  a  weak  effect  on  the 
enhancement  factor. 

There  ha  ve  been  two  st  udies  that  e  xamined  the  case  of  d  imples  much  sm  aller  th  an  th  e  height  of  th  e 
channel.  Afanaseyev6  (1993)  examined  the  case  of  a  staggered  pattern  of  shallow  CSD  with  h/d=0.067  and  /H). 25- 
0.70  in  a  channel  with  height  of  H/d=l  1-18.  Heat  transfer  was  increased  by  30-40%  with  no  measureable  increase 
in  friction  factor.  Bunker  and  Donnellan7  (2003)  examined  six  in-line  patterns  in  a  cylindrical  tube  with  h/d=0.23- 
0.4,  f=0. 34-0. 70,  and  H/d  of  approximately  4.  B  linker  reported  that  the  enhancement  was  p  roportional  to  f0'5  and 
reached  a  value  of  2  for  h/d>0.3  and  f>0.5.  The  overall  dependence  of  the  Nusselt  number  was  Re'86/"0'5. 

In  summary,  the  studies  to  date  have  found  that  the  major  variables  affecting  heat  transfer  enhancement  of  CSD 
are  h/d  and  f  and  Reynolds  number  effects  are  small.  One  stud  y  examined  the  effect  of  turbulence  intensity  and 
found  it  was  weak.  Based  on  these  parameters,  the  data  can  be  grouped  into  three  categories,  and  there  is  relatively 
little  overla p  between  them.  The  works  associated  with  Ligrani  (Moon  (2000),  Mahmood  (2002)  and  B  urgess 
(2003))  were  all  performed  at  /=0.57  and  relatively  small  values  of  channel  height,  H/d=0.2-2.  Sm  all  values  for 
channel  height  result  in  flow  di  sturbances  from  the  dimples  perturbing  the  flow  across  the  entire  channel .  Thi  s 
condition  is  relevant  to  gas  turbine  blades  but  is  not  likely  to  be  relevant  to  hot-gas  wall  heat  transfer  in  liquid  rocket 
engines.  The  studies  of  Afanaseyev  and  Bunker  examined  larger  values  of  channel  height  and  similar  ranges  of  f 
but  different  regimes  of  h/ d,  so  it  is  not  possible  to  make  direct  comparisons  between  the  data  sets.  Furthermore, 
the  correlations  that  have  been  proposed  by  Burgess  and  Bunker  have  not  included  data  from  the  other  studies  and 
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have  been  based  on  different  and  very  limited  ranges  of  f  and  h/d.  Based  on  the  literature  review  and  in  order  to 
support  the  use  of  CSD  in  expander  cycle  combustion  chambers,  the  following  research  needs  were  identified. 

1 .  The  ranges  of f  and  h/d  need  to  cover  the  entire  design  space. 

2.  Comparisons  with  previous  studies  must  be  performed  to  determine  if  the  results  are  consistent. 

3.  The  interactions  between / and  h/d  need  to  be  characterized.  A  complete  response  surface  for  heat 
transfer  is  needed. 

4.  A  study  should  be  performed  under  rocket  chamber  conditions,  because  flow  contained  within  the 
dimples  will  have  properties  significantly  different  from  the  freestream. 

5.  Large  values  of  H/d  should  be  examined. 

Based  on  these  considerations,  the  experimental  design  given  in  Table  I  was  formulated.  All  of  the  patterns  used 
a  staggered  arrangement.  The  last  column  in  the  table  is  the  increase  in  surface  area  due  to  the  addition  of  CSD.  A 
plot  of  th  e  design  points  alon  g  with  th  e  points  used  i  n  previous  stud  ies  is  g  iven  in  figure  1 .  In  ad  dition  to  th  e 
dimpled  surfaces,  a  smooth  wall  was  also  tested.  The  roughness  height  of  the  smooth  surface  was  19  microinches 
and  was  measured  using  a  diamond-stylus  profilometer. 


Table  1  Dimensional  data  for  CSD  patterns  used  in  this  study 


d  (in.) 

h  (in.) 

h/d 

/(%) 

H/d 

+Area 

#1 

0.212 

0.024 

0.111 

15 

4.7 

1.007 

#2 

0.212 

0.024 

0.111 

40 

4.7 

1.020 

#3 

0.212 

0.024 

0.111 

70 

4.7 

1.035 

#4 

0.186 

0.041 

0.222 

15 

5.4 

1.030 

#5 

0.186 

0.041 

0.222 

40 

5.4 

1.079 

#6 

0.173 

0.058 

0.333 

15 

5.8 

1.066 

#7 

0.173 

0.058 

0.333 

70 

5.8 

1.309 
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Dimple  Area  Fraction  (%) 


Figure  2  CSD  parameters  used  in  this  study  and  previous  studies 


II.  Experimental  Design  and  Methods 

Tests  were  conducted  in  the  AFRL  Aerophysics  Branch  Heat  Transfer  Rig  which  has  been  described  previously 
in  Coy  (2007)8.  For  these  tests  the  rig  consisted  of  a  gas-generator  section  with  a  2”x2”  cross-section,  a  contraction 
to  a  1  ”xl”  cross-section,  a  test  article  co  ntaining  the  test  specimens,  a  squ are-to-round  transition  and  then  a  nozzle 
with  a  throat  diameter  of  0.650”.  The  injector  consisted  of  25  shear  coaxial  elements  arranged  in  a  5x5  square  array. 
The  propellants  were  gaseous  oxygen  and  gaseous  hydrogen  and  all  tests  were  conducted  at  a  nominal  oxygen-to- 
hydrogen  mixture  ratio  of  5.5.  Each  specimen  was  tested  over  a  range  of  chamber  pressures  from  100-500  psig  in 
100  psi  increments.  A  typical  CSD  specim  en  is  shown  in  Fig.  (3).  The  test  article  is  a  heat-sink  design  with  thick- 
walled  sections  of  oxygen-free,  high-conductivity  copper.  Run  times  ranged  from  0.4  to  3  seconds.  The  rig  reaches 
quasi-steady  conditions  in  approximately  0.2  seconds. 

Fleat  flux  and  surface  temperature  were  measured  using  a  technique  based  on  measurements  of  temperature 
and  rates  of  change  of  temperature  using  embedded  thermocouples.  A  complete  discussion  of  the  theory  is  given  in 
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Coy9  (2008).  The  implementation  of  the  technique  required  designing  in  features  which  ensured  the  assumptions  of 
the  theory  were  valid.  The  thermocouples  had  to  be  installed  in  a  way  that  did  not  perturb  the  flow  of  heat,  and  the 
specimens  were  designe  d  to  ensure  the  flow  of  heat  was  one-dimensional.  T  he  test  specimens  were  fabricated  i  n 
two  halves  with  the  split  along  the  axial  centerline  of  the  test  articles.  Grooves  were  milled  into  one  of  the  mating 
surfaces  to  ace  ept  the  th  ermocouple  wires  at  a  precise  distance  from  the  hot -gas  surface.  The  hal  ves  were  braze  d 
together  with  Silvalloyl5  in  a  hydrogen  atmosphere  furnace.  Since  the  braze  joint  is  on  a  plane  of  symmetry  and  is 
approximately  0.005”  wide  and  is  itself  comprised  of  a  relatively  high  thermal  conductivity  material,  it  has  minimal 
effect  on  the  flow  of  heat.  T  o  reduce  the  e  ffects  of  a  xial  conduction,  gaps  0.010”  wide  were  included  in  the  heat 
flux  blocks  at  1”  intervals.  The  effects  of  transverse  conduction  across  the  channel  were  minimized  by  the  gaps  and 
contact  resistance  between  the  test  artic  le  and  the  specimens.  In  e  very  test  the  first  measurement  station  exhibited 
much  faster  te mperature  rise  than  th  e  downstream  lo  cations.  This  was  attributed  to  a  small  mismatch  between  the 
test  coupon  and  the  test  article  trip  ping  the  boundary  layer  and  locally  enhancing  the  heat  transfer.  Therefore  the 
data  from  the  ends  of  the  s  pecimens  were  discarded  and  only  the  center  four  transducers  were  used  to  calculate  the 
heat  flux.  Measurements  using  a  smooth  test  specimen  indicated  that  the  entrance  length  effect  in  th  e  test  rig  was 
approximately  6  inches.  This  entrance  length  corresponds  well  to  that  found  previously  by  Bunker.  The  additional 
turbulence  created  by  CSD  is  expected  to  furt  her  decrease  the  entrance  length.  However,  to  ensure  the  flow  was 
fully  developed  by  the  center  of  the  test  specimens,  an  additional  2-  inch  long  section  was  added  directly  upstream 
of  the  test  article  to  bring  the  number  of  channel  heights  at  the  measurement  location  to  6.5.  Transient  heat  flux  data 
from  a  typical  test  are  shown  in  Fig.  4.  For  this  test  the  data  used  to  obtain  the  heat  transfer  enhancement  factor  was 
obtained  by  averaging  heat  flux  over  a  data  window  from  1.8  to  2.1  seconds,  which  falls  within  the  interval  of  steady 
chamber  pressure. 

The  heat  flux  to  the  wall  depends  on  the  heat  transfer  coefficient  and  the  temperature  difference  between 
the  wall  and  bulk  fluid.  The  bulk  fluid  temperature  is  in  the  range  of  5660°  F,  so  the  percentage  change  in  heat  flux 
due  to  the  change  in  wall  temperature  is  quite  small,  approximately  5%,  so  the  heat  flux  is  quasi-steady.  The  total 
amount  o  f  hea  1 1  ransfer  t  hat  occu  rred  i  n  t  he  t  est  art  icle  up  t  o  t  he  t  hroat  can  be  est  imated  fr  om  a  C*  effi  ciency 
calculation.  For  all  o  f  the  tests  in  this  study  C*  effici ency  was  a pproximately  90%;  therefore  gas  conditions  were 
assumed  to  not  depend  on  the  CSD  pattern. 
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— HF1  — HF2  — HF3  HF4  — HF5  — HF6 


Figure  4  Heat  flux  traces,  CSD  pattern  #2,  sensor  #1  exhibits  anomalous 
heat  flux  due  to  proximity  to  leading  edge 


The  heat  fl  ux  data  fo  r  th  e  C  SD  sp  ecimens  as  well  as  for  th  e  smooth  wall  were  co  rrelated  with  Pc° 8.  This 
relationship  has  been  observed  many  times  by  others  and  has  been  codified  in  the  well-known  correlation  of  Bartz10. 
All  of  th  e  data  follows  the  expected  correlation  quite  well  except  for  the  #7  specimen  with  /=0.7  and  h/d=0.333. 
This  case  was  fit  more  closely  by  Pc° 7;  however,  for  consistency  with  the  other  cases,  the  results  presented  here  also 
assumed  Pc0'8.  Figure  5  contains  a  plot  for  several  of  the  cases  along  with  the  regression  lines.  The  R2  statistics  for 
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the  lin  ear  regressio  n  lin  es  fo  r  all  o  f  th  e  sp  ecimens  fell  bet  ween  0.992  an  d  0.  999.  The  ove  rail  heat  t  ransfer 
enhancement  factors  were  obtained  from  the  ratio  of  the  slopes  of  the  regression  lines  for  the  CSD  specimens  with 
the  sm  ooth  s  pecimen.  The  fact  that  CSD  speci  mens  ex  hibited  t  he  sa  me  Pc 0  8  scaling  as  the  sm  ooth  s  pecimen 
confirmed  the  observation  made  in  previous  studies  that  the  heat  transfer  enhancement  effect  is  nearly  independent 
of  the  Reynold’s  number.  The  results  are  given  in  Table  2  and  comprise  the  major  results  of  this  study. 


Table  2  Heat  transfer  enhancement  factors. 
Current  study. 


Specimen 

/(%) 

h/d 

Enhancement 

Smooth 

0 

0 

1.00 

#1 

15 

0.111 

1.15 

#2 

40 

0.111 

1.40 

#3 

70 

0.111 

1.48 

#4 

15 

0.222 

1.34 

#5 

40 

0.222 

1.54 

#6 

15 

0.333 

1.31 

#7 

70 

0.333 

1.82 
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•  #3 +  #4 °  smooth #1 ♦  #7 


Figure  5  Scaling  of  heat  transfer  with  Pc 8 


III.  Results 

Data  from  this  study,  as  well  as  the  two  previous  studies  with  large  H/d,  are  presented  in  Fig.  6.  Contour  levels 
have  been  fit  using  Delauney  triangulation  with  a  linear  function.  The  surface  passes  through  each  data  point  and 
uses  linear  triangulation  to  interpolate  between  the  points.  The  agreem  ent  with  the  previous  work  of  Bunker  and 
Afanaseyev  is  excellent  with  one  exception.  That  one  is  the  point  of  Bunker  with  J=0. 339  and  li  d — 0.23 3 .  This  point 
also  appears  to  not  follow  the  trend  of  the  other  data  in  Fig.  (9)  of  Bunker  (2003).  Although  we  can  cite  no  technical 
basis  for  excluding  this  point,  it  does  not  appear  to  fit  with  the  other  data,  and  therefore  we  have  chosen  to  treat  it  as 
an  outlier  and  have  not  included  it  in  the  figure.  The  fact  that  the  data  of  this  study  agrees  well  with  the  results  of 
the  earlier,  n  early-ambient  studies,  indicates  that  heat  transfe  r  enhancement  is  not  signi  ficantly  affected  by  strong 
property  gradients. 

We  can  make  a  few  observations  about  the  disposition  of  the  data  in  Fig.  6.  The  contours  at  low  values  off  are 
nearly  v  ertical,  in  dicating  that  th  e  enhancement  factor  i  s  a  fu  nction  off  only  and  nearly  independent  o  f  the  h/d 
parameter,  while  at  high  values  off  the  contours  are  horizontal,  indicating  that  enhancement  is  a  function  of  h/d  only 
and  independent  off.  In  the  mid-range  of  the  parameters  there  is  a  significant  interaction/ and  h/d. 
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Bunker  and  Donnellan  (2003)  observed  that  there  was  relatively  little  effect  off  above  a  value  of  0.5,  and 
this  ha  s  been  confirmed  he  re.  However,  the  observation  t  hat  en  hancement  i  s  p  roportional  t  o  f  '5  has  n  ot  been 
confirmed  and  may  only  be  valid  in  the  range  of  parameters  used  in  their  study. 

The  data  from  studies  with  small  values  of  H/d  are  not  included  in  Fig.  (6).  The  enhancement  factors  measured 
in  those  studies  were  in  the  range  of  2-2.5,  approximately  35%  larger  than  those  measured  here.  The  specific  trend 
with  respect  to  h/d  that  was  reported  by  Burgess  (2003)  and  given  above  as  Eq.  (1)  appears  likely  to  be  valid  only 
for  the  particular  value  off  and  H/d  used  in  that  study.  Moon’s  (2000)  result  regarding  the  effect  of  channel  height 
on  enh  ancement  is  also  not  valid  for  th  e  conditions  of  t  his  stu  dy.  Moo  n  ob  tained  the  resu  Its  fo  r  ch  annels  with 
0.37<H/d<1.49  which  may  lie  in  a  different  regime  than  H/d>4.5  considered  here. 

A  final  observation  is  that  the  increase  in  heat  transfer  is  significantly  larger  than  the  increase  in  surface  area,  as 
can  be  seen  by  comparing  Table  2  with  the  last  column  in  Table  1. 
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IV.  Summary 


This  study  measured  the  effectiveness  of  concave  spherical  dimples  for  enhancing  heat  transfer  under  conditions 
relevant  to  high-performing,  expander-cycle,  liquid  rocket  engines.  Sev  en  dimple  patterns  were  tested  and  were 
compared  tor  esults  f  or  a  s  mooth  speci  men.  The  di  mensions  of  t  he  patterns  we  re  chosen  t  o  s  pan  t  he  ra  nge  of 
possible  values  of  area  fraction  and  depth-to-diameter  ratio.  Some  conclusions  obtained  by  other  researchers  under 
near-ambient  conditions  were  confirmed.  The  enhancement  increased  with /for  all  values  of  h/d  up  to  a  value  of 
approximately  f=0. 5  and  then  did  not  increase  further,  and  the  enhancement  was  independent  of  Reynolds  number. 
However  this  study  did  not  confirm  previously  reported  measurements  of  the  scaling  of  heat  transfer  enhancement 
with  h/d  and  /  F  or  values  of  /less  than  0.2,  the  e  nhancement  became  independent  of  h/d  for  va  lues  of  h/d>  0.15; 
while  for  values  of/>0.5,  the  enhancement  continued  to  increase  with  h/d  up  to  the  highest  values  that  have  been 
tested  (h/d=0.394).  Also,  previous  studies  treated  h/d  and / as  independent  variables  while  this  study  found  a  strong 
interaction  b  etween  th  ese  v  ariables.  Fin  ally,  th  e  resu  Its  of  th  is  study  agreed  well  with  previous  stud  ies  th  at 
examined  the  case  of  channel  heights  much  larger  than  the  dimple  diameter,  H/d>4.5,  but  did  not  agree  with  a  study 
that  examined  0.37<H/d<1.49  and  found  no  effect  of  H/d;  and  this  indicates  that  there  can  be  an  effect  of  H/d  on 
heat  transfer  enhancement. 
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